Brownian colloidal particles: Ito, Stratonovich or a different stochastic interpretation 
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Recent experiments on Brownian colloidal particles have been studied theoretically in terms of 
overdamped Langevin equations with multiplicative white noise using an unconventional stochastic 
interpretation. Complementary numerical simulations of the same system are well described using 
the conventional Stratonovich interpretation. Here we address this dichotomy from a more generic 
starting point; the underdamped Langevin equation and its corresponding Fokker-Planck equation. 



Introduction. 

In a recent controversy it has been claimed, based on a 
well controlled experiment, that colloidal particles obey 
an overdamped Langevin equation with multiplicative 
noise that can not be interpreted using the standard Ito 
or the standard Stratonovich interpretations, but that in- 
stead requires an altogether different interpretation 
This particular choice has in turn been criticized and, 
based on numerical simulations, the claim has been made 
that the proposed overdamped Langevin equation must 
instead be corrected by including a new drift term, and 
that this augmented equation requires the Stratonovich 
interpretation Q . Here we will show that both positions 
are correct, but that the latter procedure can be associ- 
ated with a systematic analytical derivation from the full 
Langevin equation. 

In the experiment [l| electrically charged dielectric col- 
loidal spheres move in a solvent. The interaction forces 
have been demonstrated to be described by a conser- 
vative potential V{z) in the vertical direction z which 
includes a repulsive potential from a surface at 2; = 
and the appropriate gravitational contributions. The 
particles are also subject to thermal fluctuations and a 
position-dependent friction \{z). Both V{z) and A(z) 
have been checked and are well controlled. The particle's 
trajectory is analyzed using different refined techniques 
at nm and ms scales. 

As this is an experimentally well prepared and theoret- 
ically well described system in statistical physics, we can 
unequivocally state the main facts: the particle motion is 
driven by the classical dynamical equation of a damped 
Brownian particle and its steady state is governed by 
the Boltzmann distribution corresponding to the system 
Hamiltonian. Thus the dynamical equation is. 



■i^^\{z)z-V'{z)+g{z)m, 



(1) 



where ^(t) is a Gaussian white noise of zero mean and 
correlation function 



(2) 



obey the fluctuation dissipation relation, 

g{zf = kBT\{z)^\^{z)D^{z), 



(3) 



where the notation in [l| has been used in the last term 
to stress the relation with the diffusion coefficient D^{z). 

One can check that the stationary distribution corre- 
sponding to the Fokker-Planck equation for the proba- 
bility density Pst{z,z) is the Boltzmann distribution. 



Pst{z,z) ^ exp - 



P l2 + V{z) 



(4) 



and for simplicity we have set the particle mass to unity. 
The multiplicative function g{z) and the friction A(z) 



This is the correct starting point to analyze the ex- 
perimental observations and to clarify the theoretical 
stochastic description. 

When the friction dominates over the inertial forces 
one can approximate the second order differential equa- 
tion ([1]) by one of first order (overdamped Langevin equa- 
tion) where the acceleration term does not appear. How- 
ever, a first order Langevin equation with multiplicative 
white noise requires a stochastic interpretation [3J charac- 
terized by the value of the parameter a (the role of this 
parameter a is detailed in the appendix). The contro- 
versy occurs because of the claim in [l| that this physical 
system can be described by overdamped Langevin dy- 
namics with multiplicative noise that does not adhere 
to the well known stochastic interpretations, that is, to 
either the Ito (a — 0) or the Stratonovich (a = 1/2) in- 
terpretation. Instead, they claim that an unconventional 
stochastic interpretation with a = 1 is required. On the 
other hand, it is claimed in Q on the basis of numerical 
simulations of the same system that the correct interpre- 
tation instead requires a modification of the overdamped 
Langevin equation, and that the resulting equation must 
be interpreted according to Stratonovich, as commonly 
accepted in the physics community. 

In this brief report it will be proved that both posi- 
tions represent the same physical content, but that the 
Stratonovich proposal of '2'| can in fact be derived from 
a careful adiabatic elimination of the velocity in the full 
Langevin equation that includes the inertial contribution. 
The explicit derivation of this result makes use of calcu- 
lations concerning stochastic calculus that were carried 
out three decades ago Q. 

We next present a unified analysis of both approaches. 
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Technical details are explained in the appendix. 
The stochastic analysis. 

First it is worth to clarify that there is neither any 
doubt nor a controversy about the interpretation of the 
infradamped stochastic differential equation ([T]). We can 
operate it using the standard rules of calculus. 

The controversy starts when one wants to write the 
overdamped Langevin equation corresponding to ([T])-©. 
In this regard, it must be stressed that the overdamped 
limit is not obtained by simply setting the acceleration 
equal to zero, 



X{z) \{z) 



(5) 



which is precisely the starting point in the analysis of 
Ref. m. 

For a better comparison of the two proposals they 
can be cast in the single overdamped equation. 




V'jz) iSkeT A'(z) 



(6) 



where the correlation function of the noise is given by 
Here /3 is an auxiliary parameter to facilitate com- 
parison of the different Langevin equations under consid- 
eration. Use of ([3]) has been made to simplify the nota- 
tion. The Langevin equation proposed in [1| corresponds 
to the choice /3y = 0, while that of [H corresponds to the 
choice I3mm = 1- This is not enough and the stochastic 
interpretation of ([5]) or ([6|) matters here. 

Nevertheless the Fokker-Planck equation correspond- 
ing to a Langevin equation is a much better reference 
frame to discuss these matters. This equation has the 
following standard form, 

dP d 1 

"-Ki{z)P+-^K2{z)P, (7) 



dt 



dz 



2dz 



where Ki{z),K2{z) are the first two differential mo- 
ments, which can be calculated from the explicit algo- 
rithm ([So)) with the results (see the appendix for details) , 



K,iz) = - 
K^iz) = 2 



V'{z) ksT X'{z) 
-^-^(/3 + 2a)^; 

X{z)- 



(8) 



Then the associated Fokker-Planck equation corre- 
sponding to ^ can be written in the form. 



dP 
'dt 



d 1 



dz \{z) 

d keT \'{z) 
^d^ 2 A2(z) 



V'{z) + kBT 



d 



P 



dz 

{/3 + 2a- 2) P. 



(9) 



Now it is time to discuss both interpretations. In Ref. 
[H it is claimed that the parameter ay = 1 with /3 = 
is the correct choice, while in Ref. the choice is the 
standard Stratonovich one, ajv/M = 1/2 but with 13 — 1. 

Now it is quite clear that the two options discussed 
above cause the last term term to vanish because all of 
them fulfill the condition 



/3 + 2a = 2, 



(10) 



and accordingly they obey the same Fokker-Planck equa- 
tion. 



dP _ d 1 
'dt ~ 'd^Mzj 



V'{z) + kBT^ 
dz 



P. 



(11) 



Furthermore, this means that they reflect the same 
physics contained in the now reduced Boltzmann distri- 
bution. 



Pst{z) 



(12) 



which is a solution of pT|) . 

We can conclude that all the possible models and 
stochastic interpretations satisfying (|10l) will give the 
same physics. Moreover for those readers attached to 
the Ito interpretation, the choice would he P — 2 and 
a = 0, and the readers interested in having their own 
choice they can select a pair of /3 and a fulfilling (|10p . 

This result closes the immediate controversy, but a 
question remains open: Does a correct choice exist? The 
answer is yes, and it is based on the adiabatic elimination 
of the velocity from the full Langevin equation descrip- 
tion (HI)-©. 

The adiabatic elimination of the velocity variable 
in the full inertial Langevin equation was carried out 
following rigorous mathematical procedures in Q. The 
key point is an expansion in the A(z) as a large quantity, 
which leads to the overdamped approximation to the 
problem. This is manifest in the final form of the 
Fokker-Planck equation, which is of order \{z)~^ in the 
drift and diffusion terms. The final result is unequivo- 
cally that the physically correct overdamped equation is 
(|6]) which corresponds to the proposal of Ref. [2], that is, 
/3 = 1 along with the choice a = 1/2 corresponding with 
the Stratonovich interpretation. The technical details 
are indeed cumbersome but available to the interested 
reader, in [ij or upon request to the author. The way 
in which the final overdamped equation depends on 
the physical constrains during the adiabatic elimination 
procedure is also studied in 6]. A thermodynamic view 
of the different stochastic interpretations is discussed in 

0- 

Force and drift. 

Another discussion in P, Q is devoted to the correct 
evaluation of the force felt by the colloidal particles. If 
we know all the terms in the starting Langevin equation, 
it is obvious that the only physical force of the problem 
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comes from the potential energy f{z) = —V'{z), which 
can be extracted experimentally from the steady state 
Boltzmann distribution. Another question is the evalua- 
tion of the drift on the particle position z (which should 
not be confused with the mechanical force) in the over- 
damped description. According to most of the literature, 
the drift is just the quantity Ki{z) which appears in the 
Fokker-Planck description. In our case this is given ex- 
plicitly by 



A(z) 



A2(z) 



which depends neither on /3 nor on the stochastic inter- 
pretation parameter a because of (ITO|) . The drift can 
be obtained from experimental data on z{t) trajectories 
using its definition (see (j2ip in the appendix). 

We thus conclude that matters have been clarified 
without the need to introduce new stochastic interpreta- 
tions to understand this well controlled physical system 
either experimentally 1] or numerically [2]. The keystone 
is the connection between the infradamped equation ([T]) 
and the corresponding overdamped equation (jG)), and not 
the more intuitive overdamped equation ([S]). 
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Appendix. Here we present the derivation of the 
Fokker-Planck equation for the generic Langevin equa- 
tion with multiplicative noise [H, 



z^A{z)+B{z)at), 



(14) 



because it does depend on a parameter a e [0, 1] which 
specifies the stochastic interpretation of the Langevin 
equation. AW{t) is the Wiener increment, 



where the noise correlation is given by ©. We start 
by deriving a first order algorithm for a time step At. 
Formal integration of (IT4t leads to 

pt+At pt+At 

z{t+At) = z{t)+ / A{z{t'))dt'+ / B{z{t'))^{t')dt' 

(15) 

An expansion of the first integral gives A{z{t))At, but 
the second integral is not unequivocally defined 



t+At 



B{z{t'))i{t')dt' 



B [(1 - a)z{t) + az{t + At)] AW{t). 



(16) 



AW{t) 



t+At 



(17) 



which is a Gaussian precess with zero mean value and 
second moment 

{AWitf) = 2At. (18) 
The next step is to write an algorithm for (ITS]) to first 
order in At, which implies an expansion of the stochastic 
integral. Using the fact that to lowest order z{t + At) in 
(1161) can be approximated by 



z{t + At) = z{t) + B{z{t))AW{t) + 



(19) 



we get the explicit form of the algorithm for any stochas- 
tic interpretation a, 

z{t + At) = z{t) + A{z{t))At + B{z{t))AW{t) + 

aB{z{t))B\{z{t))AW{t)^ + .... (20) 



The Fokker-Planck equation (O for P{z,t) can be 
obtained from the Kramers-Moyal expansion [3|, where 
the quantities Ki{z) and K2{z) are obtained as follows. 
Defining the increment Az = z(t + At) — z{t), we find 



Ki(z) = lim = A(z) + 2aB(z)B'(z), 

^ ^ At^O At V ; V 

i^2(z) = lim iiMll ^ 2B{z)\ 



(21) 



These are the two quantities which appear in the Fokker- 
Planck equation ([7]). Using now the equation (jS]) we get 
the explicit expressions for A{z) and B{z), 



which are substituted in the previous equations to obtain 
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